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On the Oscillation of Combustion of a Laminar Spray 
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A spray combustor, with flow velocities in the laminar range, exhibits a unique operating mode where large 
amplitude, self-induced oscillations of the flame shape occur. The phenomenon, not previously encountered, 
only occurs when fuel is supplied in the form of fine liquid droplets and does not occur when fuel is supplied 
in gaseous form. Several flow mechanisms are coupled in such a fashion as to trigger and maintain the 
oscillatory motion of the flame. These mechanisms include heat transfer and evaporation processes, dynamics 
of two-phase flows, and effects of gravity (buoyancy forces). An interface volume, lying between the fuel 
nozzle and the flame was found to be the most susceptible to gravity effects, and postulated to be responsible 
for inducing the oscillatory motion. Heptane fuel was used in the majority of the tests. ^ ^ ^ '7 


INTRODUCTION 

The classical paper of Burke and Schumann [1] 
on gaseous diffusion flames is considered a 
milestone in the analysis of combustion pro- 
cesses. Although the authors made a number 
of simplifying assumptions in their model, in- 
cluding the neglect of gravity, constant gas and 
air velocities, and a constant diffusion coeffi- 
cient for all species, their predictions were 
found to show relatively good agreement with 
experimental measurements. Gas composition 
and flame dimensions were predicted at vari- 
ous locations above the fuel nozzle. Other re- 
searchers [2-4] followed in their footsteps and 
improved the original model, but always re- 
tained the assumption that the fuel was sup- 
plied in gaseous form and neglected gravita- 
tional effects. 

Various observations of a Burke-Schu- 
mann-type flame under laminar operating con- 
ditions revealed the phenomenon of flame tip 
flickering at frequencies of about 10-20 Hz. 
This effect, which was observed experimentally 
[5-9] and also predicted numerically [10, 11], is 
the direct result of gravity-induced buoyant 
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forces on the flow of the reaction products at 
some distance downstream of the nozzle, aero- 
dynamic instabilities in the jet itself and heat 
release effects. Ellzey et al. [11] presented pre- 
dictions both considering and neglecting grav- 
ity. They showed that neglecting gravity gave a 
steady flowfield but including gravity caused 
large structures to form that were convected 
downstream at frequencies on the order of 
15-20 Hz. An inviscid stability analysis for a 
buoyancy-driven flow in an infinite candle [12] 
predicted the frequency of the flickering within 
a factor of 2. Hence, it appears that the rele- 
vance of gravity is significant for these types of 
flames. 

Some work has been previously reported for 
laminar flames using fuel supplied as a liquid. 
A theoretical and experimental study [13] of a 
laminar diffusion spray flame did not en- 
counter instabilities in the flame. Gravity was 
considered in this model, however, the model 
was derived only for steady-state conditions 
and consequently no instabilities were pre- 
dicted. In addition, the liquid fuel was supplied 
directly from a nozzle as a conical spray that 
was surrounded by the combustion air, which is 
different from the Burke-Schumann configu- 
ration. 

Studies by Greenberg [14-16] investigated 
the effects of fuel supplied as a liquid on the 
behavior of laminar diffusion flames in the 
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Burke-Schumann configuration. In these stud- 
ies, the central gaseous fuel jet used in the 
original formulation was replaced by liquid fuel 
droplets of various diameters carried by a flow 
of nitrogen. Effects of various initial velocity 
distributions, injection angles for the droplets, 
and preevaporation of the liquid fuel droplets 
were investigated. However, gravity was not 
included in these studies, and no instabilities 
in the flame behavior were predicted. An un- 
steady behavior of the flame was predicted 
[17], when gravity was neglected in the formu- 
lation, only when the fuel was supplied in a 
cosinusoidally alternating fashion of liquid and 
gas. 

The present experimental study was moti- 
vated by the analytical study of Greenberg [14, 
15] and Greenberg and Cohen [16, 17]. The 
objective of the present study was to experi- 
mentally investigate the effects of fuel droplet 
diameter, preevaporation of the liquid fuel and 
droplet trajectories on the structure of the 
flame in a Burke-Schumann configuration. The 
combustor utilized for the present study was a 
coaxial flow system flowing in the upward di- 
rection. The inner stream (fuel stream) con- 
tains fuel in the firm of liquid droplets and fuel 
vapor in different proportions. The fuel was 
carried in the upward direction by a nitrogen 
gas stream. Air flowed in the outer coflowing 
stream. Both the fuel and air streams had 
equal average velocities. 

EXPERIMENTAL APPARATUS 

The combustor was designed to study laminar 
spray combustion in an extension of the classi- 
cal Burke-Schumann diffusion flame configu- 
ration. A schematic drawing of the combustor 
is illustrated in Fig. 1. The combustor consisted 
of two coaxial tubes. The outer tube was quartz 
and had an i.d. of 71 mm. Two different inner 
tubes were utilized in the present study. A 
quartz inner tube was used for internal laser 
diagnostics, and had an i.d. of 12.90 mm and an 

0. d. of 14.80 mm. A stainless steel tube with an 

1. d. of 12.54 mm and an o.d. of 15.80 mm was 
also used in the study when laser diagnostics 
were not required. The length (height) of the 
outer tube was 30 cm. Tests with 100-cm and 
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Fig. 1. Schematic drawing of the combustor. 


10-cm length tubes or even without the outer 
tube, were performed to eliminate the possibil- 
ity of acoustic coupling effects and showed 
negligible differences in behavior. As illus- 
trated in Fig. 1, the combustion air entered 
through four radial entrance holes and passed 
through a series of fine-mesh screens and flow 
straighteners. The upstream length from the 
screens to the combustor (the tip of the center 
tube) was about 25 cm and the total length of 
the center tube (from the atomizer to the tip) 
was about 50 cm. 

Fuel was introduced into the combustor 
through the center tube in the form of fine 
droplets. A slightly modified Sonotek ultra- 
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sonic nozzle #8700-60 MS or a self-assembled 
ultrasonic unit were used to generate the 
droplets. As shown in Fig. 2, the droplets were 
injected upwards as a narrow conical spray. 
The droplets were then carried upwards by a 
nitrogen gas stream. The nitrogen was injected 
through a tube which surrounds the “horn” of 
the ultrasonic atomizer. As the droplets trav- 
eled upwards through the 50-cm tube, they 
partially vaporized and interacted with each 
other, causing a subsequent change in their 
size and momentum. The large droplets were 
carried upwards at a velocity proportional to 
the difference between the nitrogen carrier gas 
flow and their Stokes velocity. The small 
droplets (less than 4 microns) can be used to 
represent the carrier gas velocity. This is illus- 
trated in Fig. 3. Figure 3 presents centerline 
phase/Doppler data taken at axial locations of 
300 and 20 mm below the exit of the center 
tube. Results are presented for droplet velocity 
and normalized drop number for droplet diam- 
eters up to 60 fim. Experimental uncertainty 



in measurements for drop size is estimated at 
±6.5% and ±0.5% for velocity. Flowrates of 
the nitrogen and fuel streams were maintained 
at 13.2 X 10 “ 6 and 12.0 X 10 6 kg/s, respec- 
tively. Time resolved measurements of both 
the gas and drop velocities for up to 10 s inside 
the center tube indicated both to be steady. 
Heptane fuel was used for the results reported 
in the present paper. 

Figure 4 presents velocity profiles for three 
drop sizes across the center tube at 20 mm 
upstream of the center tube exit at nearly the 
same conditions as Fig. 3. Flowrates of the 
nitrogen and fuel streams were 13.2 X 10' 6 
and 14.0 X 10 -6 kg/s, respectively. The rela- 
tively low velocity of the carrier gas near the 
wall of the inner tube was not large enough to 
carry the larger drops. This caused them to 
descend and come into contact with the tube 
wall forming a liquid film that flowed down- 
ward. The number of droplets that hit the wall 
and the resulting flow in the downward direc- 
tion was significant, especially in the lower part 
of the center tube. A drain was incorporated in 
the combustor configuration at the bottom of 
the center tube for continuous removal of the 
liquid film. The measurement of the actual fuel 
flow rate to the combustion zone always con- 
sidered the drained amount. 

Experimental measurements were obtained 
using two techniques. These included a two- 
component phase/Doppler particle analyzer 
(PDPA), which simultaneously measured the 
size and velocity of individual droplets, and a 
CCD camera coupled with laser light-sheet 
illumination for large-scale observations of the 
two-phase flow and flame. The average flame 
oscillation frequency was measured by the vari- 
ation in gas velocity over a number of cycles in 
the combustion zone using measurements from 
the PDPA system. 

RESULTS AND DISCUSSION 

The oscillatory motion encountered in the pres- 
ent study is different from any reported previ- 
ously in the literature. In this study, the entire 
volume of the flame was affected in an expan- 
sion/contraction type of modulation, whereas 
previously reported flame flickering affected 
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Fig. 3. Variation of centerline vertical 
velocity and number density with drop 
size at two axial locations inside the 
center tube. 


the flow of reaction products only at distances 
downstream of the nozzle. 

Two sets of pictures are presented in Figs. 5 
and 6. Figure 5 illustrates the flame in a typical 
cycle of the phenomenon displayed in a com- 
plete view of the combustor and flame. The 
outer quartz tube and the inner stainless cen- 
ter tube are clearly visible. The flame is indi- 
cated by a blue color with a glowing soot tip. 
Flowrates of the air, nitrogen, and heptane 
streams were 107 X 10 b , 10.1 X 10“ 6 , and 
2.84 X 10" 6 kg/s, respectively. Pictures are 
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Fig. 4. Radial distribution of vertical velocity 20 mm up- 
stream of the center tube exit. 


marked by a video clock running at 30 Hz, 
hence each clock increment corresponds to 
33.3 ms. For example the specific flame shown 
has an amplitude ratio (average ratio of mini- 
mum to maximum flame height) of about 1:11 
and an average frequency of 3 Hz. The fre- 
quency was obtained from time-resolved PDPA 
measurements over a number of cycles. A much 
closer view of the phenomena is given in Fig. 6, 
where the camera was focused to image closer 
to the exit of the center tube containing the 
fuel droplets and nitrogen. In Fig. 6, the 
flowrates of the air, nitrogen, and heptane 
streams were 119 X 10“ 6 , 11.8 X 10” 6 , and 
13.5 X 10 -6 kg/s, respectively. It is evident 
that the equivalence ratio in Fig. 6 is much 
larger than in Fig. 5, 1.7 versus 0.4, however, 
the basic phenomena remained essentially 
identical except for absolute flame dimensions 
and frequency of oscillation. For flow visualiza- 
tion, the coflow air stream was seeded with 
aluminum oxide particles with a nominal size 
of 1 /Am. The complete flowfield was illumi- 
nated by a light sheet (from left to right in the 
pictures) from an Ar + laser. The particles ob- 
served close to the center tube tip are the fuel 
droplets since the nitrogen stream was not 
seeded. The center tube was transparent quartz 
in order to provide optical access. The high- 



Fig. 5. Complete cycle of flame phenomenon. Test conditions: 107 x 10- 6 kg/s air, 10. 1 x 10- 6 kg/s nitrogen, 
and 2.84 x 10-* kg/s heptane. 










Fig. 5 Continued 







Fig. 6. Closeup view of flame phenomenon cycle with laser-sheet illumination and coflow air stream seeded. 
Test conditions: 119x10-6 kg/s air, 1 1 .8 x 10-* kg/s nitrogen, and 13.5 x 10-6 kg/s heptane. 
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intensity illumination from the laser masks the 
blue flame color, however, the sooting tip is 
still clearly evident. As previously discussed, 
the flow conditions in Fig. 6 are not identical 
to those in Fig. 5; consequently, the average 
frequency shown here is 2.5 Hz. 

Except for two extreme cases of either very 
low velocities (up to about 0.01-0.02 m/s) or 
relatively high velocities (above about 0.3 m/s) 
of the nitrogen and air streams, the flame 
always exhibited an oscillatory motion when 
the fuel was supplied as a liquid. Both the 
flame dimensions and the frequency of the 
oscillations varied depending on the upstream 
flow conditions and ranged from 1-5 Hz for 
the frequency and 1:20 in amplitude ratio for 
the flame dimensions. 

Following many sets of such photographs 
and local instantaneous size and velocity mea- 
surements performed with the phase/Doppler 
system, it is possible to establish a preliminary 
description of the phenomenon. 

The oscillatory motion of the flame appears 
to originate from several coupled thermal and 
fluid mechanic mechanisms. A preliminary de- 
scriptive explanation of the mechanism and 
sequence of events within a cycle is proposed. 

At the start of the cycle, before combustion 
commences, the fuel stream containing nitro- 
gen, fuel droplets and fuel vapor flows upward 
in the center tube until it reaches the top of 
the center tube. After it emerges from the 
center tube, the fuel stream immediately ex- 
pands radially outward and then flows down- 
ward creating a mushroom shaped cloud of 
nitrogen, fuel droplets and vapor. This is due 
to the significant density difference between 
the fuel stream and the surrounding coflowing 
air. The fuel stream is ignited by a descending 
flame which initiates the reaction between the 
fuel vapor and the adjacent oxygen (see Figs. 
5b and 6a). After ignition, the flame rapidly 
expands predominately in the radial direction 
(see Figs. 5c, 5d, 6b, and 6c), consuming the 
fuel vapor that exists at the interface between 
the top of the fuel stream and the coflow air. 
Thereafter, the flame spreads at a slower rate 
only in the vertical direction and the rate is 
now controlled by the evaporation rate of the 
fuel droplets (see Figs. 5e-5j and 6d-6i). Buoy- 
ancy effects lift the reaction products and, at a 


different rate, the emerging fuel stream. The 
fuel stream, which originally flowed radially 
outward and downward after emerging from 
the center tube, now acquires a vertical veloc- 
ity component (see Fig. 6d-6i). This change is 
primarily due to heat transfer from the flame 
to the fuel stream by conduction and radiation. 
The height above the center tube that the fuel 
stream reaches appears to be a critical param- 
eter in the dynamics of this flame. This height 
is determined by the velocities of the coflow 
and fuel streams. 

The oxygen for reaction is available from 
two sources. The first source is radial diffusion 
inward from the surrounding coflow. This dif- 
fusion process must overcome the radial out- 
ward velocity of the newly generated reaction 
products. Oxygen is also available from a por- 
tion of the coflow air stream that manages to 
penetrate the combustion zone from the bot- 
tom of the flame, close to the center tube lip. 
This flow of oxygen is visualized by the laser- 
light scattered from the aluminum oxide seed- 
ing particles in the coflow air stream and is 
best seen in Figs. 6d-6f. 

The fuel is consumed from the outer portion 
of the fuel stream nearest the flame, causing a 
reduction in the length and radius of the fuel 
stream (Figs. 6g, 6h). This rapid reduction in 
the radial dimension induces a short swirling 
motion (Figs. 6h, 6i), probably due to a tempo- 
ral instability in the flow. The swirling motion 
is difficult to visualize from the photographs 
presented, distinguished only by horizontal 
lines from particle trajectories (Figs. 6h, 6i), 
however, it is quite evident when observing a 
videotape of the sequence of events. The con- 
sumption of fuel vapor in the fuel stream leads 
to a reduction in flame size (see Figs. 5j-51 and 
6j— 61). This reduction in the flame size reduces 
the amount of heat transferred from the flame 
to the emerging fuel stream. Consequently, the 
fuel stream temperature is reduced causing a 
decrease in vertical velocity until negative val- 
ues are observed (Figs. 6i-6k). A separation 
develops between the flame and the top of the 
decelerating fuel stream that prevents a fresh 
supply of fuel vapor from reaching the flame 
(see Figs. 6k, 61). 

The upward coflowing air now faces a de- 
scending semirigid center bluff body consisting 
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of the relatively dense, droplet-laden fuel 
stream. Consequently, a toroidal recirculation 
zone is formed in the wake of the fuel stream 
with its center moving downwards (Figs. 61, 
6m). This toroidal recirculation increases the 
radial spreading of the fuel stream (Figs. 6m, 
6n), reduces the position of the flame due to 
the reduced pressure at its center (Figs. 6a, 6b, 
61, 6m) and transfers fuel vapor from the top of 
the fuel stream upwards due to the upward 
velocity at the outer edge. The flame descends 
until it approaches the top of the emerging 
fuel stream which contains a large quantity of 
fuel vapor from vaporization of the droplets 
(Figs. 6b, 6n, 5b, 5m). The flame then expands 
rapidly again in the radial direction and the 
sequence of events is repeated from the begin- 
ning. 

Hence, it appears that the modulation in the 
flame shape originates from the coupling be- 
tween the flame and the fuel source. Although 
the fuel is supplied through the center tube in 
a steady fashion, confirmed experimentally with 
phase/Doppler measurements near the exit of 
the center tube, the region between the exit of 
the center tube and the flame is sufficient to 
allow for the fuel stream to periodically expand 
vertically and collapse radially, thus inducing 
the unique oscillatory motion observed. This 
phenomenon is altogether different from the 
flame flickering previously reported in the lit- 
erature due to the large toroidal vortices gen- 
erated by the buoyantly induced flow. 

As previously discussed, at very low veloci- 
ties, no oscillatory motion was encountered. In 
this case, the carrier gas velocity was sufficient 
to carry only the very small droplets, while the 
larger ones accumulated on the inner wall and 
drained out of the system. Because the small 
droplets evaporated almost completely prior to 
their entrance into the combustion zone, the 
fuel was essentially supplied in gaseous form 
and stable combustion was achieved. At the 
other extreme, higher velocities on the order 
of 0.3 m/s also established a quasi-steady 
flame. For this case nearly all the generated 
droplets were carried upwards to the combus- 
tion zone. The resulting flame height was rela- 
tively constant with some random oscillations 
of smaller amplitudes and at higher frequen- 
cies. At the higher flowrates, the increased 


momentum of the coflow air stream prevented 
the collapse of the fuel stream as it emerged 
from the center tube. This increased the height 
above the center tube that the fuel stream 
reached providing increased residence time for 
droplet evaporation to occur. The higher veloc- 
ities also increased the mixing of fuel vapor 
and air. In addition, at the higher velocities, a 
region of low velocity was formed in the wake 
of the central fuel stream that served to stabi- 
lize the flame. 


CONCLUSIONS 

A self-induced and continuous oscillatory mo- 
tion, observed in liquid spray combustion with 
low velocities, is described. This phenomenon 
was observed over a relatively wide range of 
operating conditions and practically excludes 
the possibility of achieving laminar conditions 
in similar configurations. 

A partially stable spray flame was achieved 
only at relatively high velocities where some 
turbulence in the flame structure was ob- 
served. A stable flame was also observed at 
extremely low velocities, but all droplets had 
vaporized by the time they reached the com- 
bustion zone. A preliminary descriptive expla- 
nation of the processes involved and support- 
ing evidence are given. 
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